Distributed cavity phase frequency shifts of the 
caesium fountain PTB-CSF2 



S. Weyers\ V. Gerginov^ N. Nemitz\ R. Li^ and K. Gibble^ 

^ Physikalisch-Technische Bundesanstalt (PTB), Bundesallee 100, 38116 
Braunschweig, Germany 

^ Department of Physics, The Pennsylvania State University, University Park, PA 
16802, USA 

E-mail: stefan.weyers@ptb.de 
Abstract. 

We evaluate the frequency error from distributed cavity phase in the caesium 
fountain clock PTB-CSF2 at the Physikalisch-Technische Bundesanstalt with a 
combination of frequency measurements and ab initio calculations. The associated 
uncertainty is 1.3 x 10~^^, with a frequency bias of 0.4 x 10~^^. The agreement between 
the measurements and calculations explains the previously observed frequency shifts 
at elevated microwave amplitude. We also evaluate the frequency bias and uncertainty 
due to the microwave lensing of the atomic wavepackets. We report a total PTB-CSF2 
systematic uncertainty of 4.1 x 10~^^. 
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1. Introduction 

Over the past decade the number of primary caesium fountain clocks contributing to 
International Atomic Time (TAI) has steadily increased and now exceeds the number of 
contributing primary caesium beam clocks. During the last two years, ten fountain 
clocks from seven laboratories contributed to TAI [Ij. The fountain clocks largely 
calibrate TAI because of their excellent systematic and statistical uncertainties. While 
total systematic uncertainties clearly below 10~^^ are routinely reported by several 
fountains, there has not been widespread agreement about the methods to establish 
the frequency uncertainties due to spatially distributed phase variations of the field in 
the microwave cavity. 

Until recently, distributed cavity phase (DCP) was the largest uncertainty 
contribution for several fountains [21 |3l HI E]. The resistive losses of the cavity walls 
lead to small traveling waves, which account for the power flow from the cavity feeds to 
the walls and produce the spatial phase variations [6]. A recent comprehensive model 
has clarified the evaluation of DCP frequency shifts [3]. The ah initio model, with no 
free parameters, has been stringently verified with the French fountain SYRTE-F02 
and already used to evaluate the DCP uncertainties of it and the British NPL-CsF2 
fountain [H|5]. Here we report the results of measurements and calculations, using this 
model to evaluate the DCP uncertainty of the PTB-CSF2 fountain clock. In contrast 
to SYRTE-F02 and NPL-CsF2, the amplitudes and phases of the two cavity feeds of 
PTB-CSF2 are balanced by symmetry and cannot be independently adjusted [3 [8]. 
Here, we also evaluate the microwave lensing frequency shift [5l [9] and summarize all of 
the evaluated systematic biases and uncertainties of PTB-CSF2. 

The fountain geometry naturally minimizes frequency shifts from phase variations 
of the microwave field. In atomic fountain clocks, the atoms interact twice with the 
microwave field in the same cavity, eliminating the end-to-end phase shifts of beam 
clocks [ini [H] with a velocity reversal, which largely cancels the Doppler shifts from 
longitudinal phase gradients. Transverse phase variations clearly produce frequency 
shifts when the atoms probe the field with different phases when they traverse the 
cavity at different positions during their ascent and descent in a fountain. 

The central feature of the model P, [3] is to decompose the transverse phase 
variations into a Fourier series cos(m0) with cylindrical coordinates {p,(j),z). This 
enables densely-meshed finite element calculations of the cavity fields because the 3D 
problem is reduced to a short series of 2D solutions in p and z. The transverse phase of 
the field is proportional to p'"cos(m0) to lowest order and therefore only m < 2 terms 
are significant since the cavity apertures restrict the atomic trajectories to small p. The 
azimuthal symmetry of each Fourier term specifies the physical effects in the fountain 
that produce DCP frequency shifts. While the m = transverse variations are small, 
the m = 1 phase variations are naturally large, representing a linear phase gradient and 
a power flow from a cavity feed across the cavity to the walls. This gradient is highly 
suppressed in all fountains by using two or more diametrically opposed feeds. The 
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quadrupolar m = 2 transverse phase variations are maximally excited by two opposing 
feeds and can be negligible with four (or more) azimuthally distributed feeds [12]. 

The DCP shifts of the lowest Fourier terms usually have unique functional 
dependences on the clock's microwave amplitude. For example, if the transverse phase 
variations have no longitudinal variation for all atomic trajectories, then each atom sees 
the same phase for the entire cavity traversal during ascent and descent. Then, the 
DCP frequency shifts are independent of microwave amplitude. But, in general, there 
are longitudinal phase variations and therefore the atoms see a time-dependent phase 
during each cavity traversal, which often gives a large microwave amplitude dependence 
of the effective phase of the cavity field |3l|13]. Even longitudinal phase gradients that 
have no transverse phase variation can produce large amplitude-dependent frequency 
shifts because the atoms do not have a transverse velocity reversal in the fountain 
and therefore experience different average microwave pulse areas on the two cavity 
passages [3l [13]. The m = 0, azimuthally symmetric phase deviations have a neghgibly 
small transverse variation and a very large longitudinal variation, primarily due to the 
endcap losses [6]. Here we measure the resulting large, m = amplitude-dependent 
frequency shift, explaining the previously observed frequency shifts of the PTB-CSFl 
and PTB-CSF2 fountains at elevated microwave amplitudes [T^t [15]. 

2. DCP evaluation of PTB-CSF2 

The cylindrical TEqh Ramsey cavity of PTB-CSF2 is made from OFHC copper, with 
an inner diameter of 48.4mm and 28.5 mm height [71 [8]. In the top and bottom endcaps 
two apertures 10 mm in diameter enable the atoms to pass through the cavity, with 
70 mm long waveguide cutoff tubes to minimize microwave leakage from the cavity. The 
coaxial line that feeds the microwave signal is strongly coupled to a standing wave 
in a curved waveguide. Two opposed slits (10.8 mm length and 2 mm width) in the 
curved waveguide feed the cavity. Their strong coupling to the cavity gives a relatively 
low loaded quality factor Q = 1600, which reduces the effect of the cavity detuning, 
for example due to temperature changes, as compared to weak coupling. Our finite- 
element calculations O [3] of the cavity geometry give consistent resonant frequencies 
for the measured TEqh and TMm modes. 

In PTB-CSF2 the entire atom collection and cooling chamber, including the laser 
beam collimators, can be pivoted around its center on a goniometric stage. It allows a 
sub-milliradian adjustment of the launch direction to maximize the number of detected 
atoms. The initial center-of-mass position of the atomic cloud is offset from the 
fountain axis by 1.7 mm, along the axis of the two cavity feeds. This is inferred from 
the single Rabi pulse spectrum of \F = 3,mp = 0) — )■ \F = 4,mF = ±1) transitions, 
taking advantage of a slight unintended tilt of the static magnetic field in the vertical 
plane perpendicular to the feeds [16]. Because of this offset, and since PTB-CSF2 is 
vertically adjusted to maximize the number of detected atoms, the actual angle a\\ at 
which the detected atoms are launched deviates from true vertical by 0.15 mrad so that 
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the atomic cloud is launched with a nonzero horizontal launch velocity of 0.7mm/s in 
opposite direction of the offset. Perpendicular to the feed axis, the above evaluation 
of the spectrum cannot be used since the static magnetic field happens to have a 
negligible component parallel to the feeds [16]. Instead we evaluate the |3,0) — )■ |4, 1) 
transition probabilities for a single Rabi pulse on the ascent as a function of the launch 
angle perpendicular to the feed axis. This transition probability shows no gradient 
at the launch angle that maximizes the number of detected atoms for a cloud that is 
launched on the cavity axis. We observe no gradient of the transition probability versus 
perpendicular tilt and therefore conclude that there is no cloud offset perpendicular 
to the feed axis. For tilts parallel to the feed axis, this technique is consistent with 
the above measured offset because there is a gradient of the |3,0) — )■ |4, 1) transition 
probability versus launch angle. 

The validity of the Monte-Carlo model of the fountain is confirmed by comparing 
experimental and calculated transition probabilities P as a function of microwave 
amplitude b (see figure [I]) [H [17]. The agreement is sensitive to the size, offset, and 
velocity distributions, the positions of apertures in the fountain, and the detection 
inhomogeneity. Here b is in units of 7r/2-pulses during a cavity traversal for a uniform 
density distribution on a cavity passage [3]. Because the cloud size in PTB-CSF2 is 
smaller than the cavity aperture during the first cavity passage, the optimum microwave 
amplitude for normal clock operation, which maximizes the Ramsey fringe contrast, is 
b = 0.97. 

We next evaluate the DCP frequency shifts from each Fourier component m = 
0,1,2. The frequency shifts are often singular near even b where the Ramsey fringe 
contrast vanishes. Therefore, we instead discuss the change in transition probability 
6P at detuning +Az//2, where Az/ = 0.912 Hz is the full width of the central Ramsey 
fringe. The frequency shift is 6P/ (dP/du) where the Ramsey fringe slope for PTB-CSF2 
is dP/dv = 1.57/Hz at optimum amplitude. 

2.1. m = DCP frequency shift 

The m = DCP shift is primarily due to power flow from the feeds at the cavity 
midplane to the endcaps. This power flow causes large longitudinal phase variations, 
which are symmetric about the cavity midplane, and small, azimuthally symmetric 
transverse phase variations. The longitudinal phase variations are not canceled by the 
velocity reversal of the two cavity passages because the cloud is larger for the descent 
than the ascent, and therefore experiences a different average fleld amplitude [31 IT3] . 
This gives large frequency shifts at elevated microwave amplitude, especially 6 = 4 and 
8 where the sensitivity function gives the maximum perturbation for a phase variation 
that is symmetric about the cavity midplane. At low amplitude, 6 <^ 1, 6P is small 
because the effective phase is zero since it is simply the average phase along a trajectory 
through the cavity. At 6 = 2, 6P is also zero because the phase is an even function 
of 2; [21 |13]. Thus, 6P is highly suppressed at optimum amplitude because, while it 
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increases from zero as for small 6, it cannot get too large since it returns smoothly to 
zero near 6 = 2. 

The solid curve in figure [2] shows the calculated changes 5P of transition probability 
as a function of microwave amplitude for PTB-CSF2. The initial measured l/\/e cloud 
radius at launch is 3.9 mm and the temperature is 1.2 /iK. If the cloud is initially centered 
(no offset), the 5P is almost identical. The curves show the characteristic features of 
m = DCP shifts: a very small shift around optimum amplitude 6 < 1 and large shifts 
near 6 = 4 and 6 = 8. At 6 = 8 the m = DCP shift is 7.9 x 10"^^ 

We extract the experimental change in transition probability in figure [2] (full circles) 
from measurements of the frequency difference between the PTB-CSF2 and PTB-CSFl 
fountains, where PTB-CSFl [18] is operated at optimum amplitude. We also measure 
the Ramsey fringe slopes of PTB-CSF2 at each b. We apply frequency corrections 
to both fountains for the quadratic Zeeman effect, cold collisions, blackbody radiation, 
gravity and the relativistic Doppler shift. The error bars are 1 a statistical uncertainties. 
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Figure 2. The to = DCP change of transition probabihty 5P in PTB-CSF2 as 
a function of microwave amphtude h. The initial cloud offset is 1.7 mm. Solid line: 
Calculated 5P for an initial cloud radius at launch of 3.9mm [1/^). Gray band: 
Range of shifts for initial cloud sizes of ±20%. Full circles: Measurements. 

At optimum amplitude of PTB-CSF2, h = 0.97, the two fountains agree within the 
statistical uncertainty. 

The measurements in figure [2] must also include m = 1 and m = 2 DCP shifts. 
Below we show that both are negligible on this scale (see sections 12.21 and [23]) . The 
calculated m = DCP shifts agree well with the measurements, which explains the 
previously observed amplitude dependence of the PTB-CSF2 fountain [7j. 

At optimum amplitude the calculated frequency shift is 1.8 x 10^^®. Because the 
m = DCP shift depends strongly on the difference of the average microwave amplitudes 
that the atoms experience on the two cavity passages, the shift is sensitive to the initial 
cloud size. In figure [2] the gray band shows the range of shifts for initial cloud sizes of 
±20%. This range of initial cloud size gives a negligible uncertainty of 6 x 10~^^. In [3], 
measured m = 1 DCP shifts suggested that the cavity walls could have inhomogeneous 
conductivities as large as 20%. If the conductivity of the top and bottom endcaps differ 
by 20%, it would give an m = DCP shift of 3.8 x 10~^® at optimum amplitude. The 
measurements in figure [2] do not exclude this possibility and we therefore take this as 
our m = DCP uncertainty. 
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2.2. m = 1 D CP frequency shift 

The m = 1 DCP shift is caused by power flow across the cavity, arising from potential 
imbalances between opposing feeds or resistance inhomogeneities of the copper cavity 
walls. Near the cavity axis, the m = 1 phase variations are linear gradients and lead 
to frequency shifts if the mean transverse positions of the atomic cloud are different on 
the two cavity passages, for example if the fountain is intentionally tilted [HI |3l IH [5] . 
Because power is supplied to the PTB-CSF2 cavity via a single cable to the curved 
waveguide, we cannot experimentally amplifiy the m = 1 DCP shift by alternately 
supplying one feed or the other when the fountain is tilted [HI HJ |5]. We therefore 
cannot use this technique to balance the feeds or to measure the tilt of the fountain to 
give no atomic motion along the feed axis and no m = 1 DCP shift. 

To evaluate the m = 1 DCP shifts, we measure the frequency change of PTB- 
CSF2 at optimum amplitude as we change the launch angle to ±2.54 mrad from a\\ 
and a± that maximize the number of detected atoms, parallel and perpendicular to 
the feed direction (see figure [3]) [HIS]. Here, we again use the PTB-CSFl fountain as 
the reference. Phase gradients parallel to the feeds can be caused by feed imbalances 
and resistance inhomogeneities, while phase gradients perpendicular to the feeds can 
only be from resistance inhomogeneities. The measured tilt sensitivities are my = 
1.9(8.7) X lO^^^mrad^^ for parallel and m± = 7.2(10.7) x 10^^''mrad^^ for perpendicular 
tilts. 

We need to determine the difference between the launch angle a\\ for normal 
operation and the launch angle an^O) where the m = 1 DCP shift is zero. Note that 
cancelling the m = 1 DCP shift of an offset cloud implies launching the atoms towards 
an even larger mean offset for the descending cloud since it in effect uniformly fills the 
cavity aperture. Therefore, a moderate offset implies that significantly fewer atoms 
are detected. For the same reason, a vertical launch gives a non-zero m = 1 DCP 
shift because the ascending cloud is offset and the descending cloud is nearly centered. 
A Monte-Carlo simulation shows that the 1.7 mm initial cloud offset corresponds to 
a launch angle difference of a\\ — ay^o = 1.09 mrad. Its uncertainty cToJi = 0.57 mrad 
includes the quadratic sum of the uncertainty of the mean position (0.27 mrad) [16] and 
possible variations (0.5 mrad) from long-term molasses laser-beam intensity imbalances. 
From m||(a||o — ay), the m = 1 DCP shift is +2.1 x 10"^'' and its uncertainty is 
H<|| + <||<|| + <||(«ll,o - «||)^)^/^ = 10.8 X 10-1^ 

Perpendicular to the feeds, the cloud offset is consistent with zero (see section [2]) 
so that a vertical launch maximizes the number of detected atoms. This gives no 
m = 1 DCP shift for atomic motion perpendicular to the feeds. The |3,0) — )■ |4, ±1) 
measurements yield an angular uncertainty of 0.33 mrad. Here, we again add in 
quadrature an 0.5 mrad uncertainty related to possible long-term molasses laser-beam 
intensity imbalances. Thus, the m = 1 DCP uncertainty for the perpendicular direction 
is {m']_a'^^^ + o"^,_L^a,±)^^^ = ^-^ ^ 10"^''. For both tilt directions, we get a total m = 1 
DCP shift of 2.1(13.3) x 10"^^. 



Distributed cavity phase frequency shifts of the caesium fountain PTB-CSF2 



in 
O 



CO 

o 

CNI 
LL 

O 



1 - 







1 1 

(a) 1 


1 


1 
















1 





-3-2-10123 




Tilt angle / mrad 



Tilt angle / mrad 



Figure 3. CSF2 frequency at optimum amplitude versus tilt, ±2.54 mrad, parallel 
(a) and perpendicular (b) to the feed direction. Zero tilt (dashed line) corresponds 
to the launch angles a|| and a_L that maximize the number of detected atoms. The 
dash-dotted line in (a) indicates the launch angle a||_o a-t which the m = 1 DCP 
contribution in the feed direction vanishes. The weighted fits (solid lines) are (a) 
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1.9(8.7) X 10-"mrad-i and (b) mj. = 7.2(10.7) x IQ-^mrad" 



It is interesting to deduce the corresponding feed imbalance from the measured 
parallel tilt sensitivity in figure [3]^a). In figure H] we show the calculated microwave 
amplitude dependence of the m = 1 DCP transition probability change 5Pw of PTB- 
CSF2 for a single weakly coupled feed. The most likely imbalance for our cavity is an 
asymmetry in the machined size of the cavity couplings. Such an asymmetry supplies 
more power to, but also leads to greater loss from one feed. Therefore, a feed imbalance 
of e gives 5P = e/2 x [3J. A 10% imbalance yields a 5P of only 5% of 5Pw in figurelH 
If, however, the microwave amplitude in the curved waveguide is imbalanced at the two 
couplings, 5P is larger by a factor Qo/Q [3], where Qo is the theoretical quality factor for 
weak coupling, calculated to be 28600 for PTB-CSF2, giving Qo/Q = 17.9. Assuming no 
resistance inhomogeneities, the tilt sensitivity in figure |3]^a) only constrains the couplings 
to be the same to 36% and the microwave amplitudes in the curved waveguide to 2%. 
We expect the balance of both is better than these estimated limits. Additionally, 
the measured m = 1 DCP shift of 2.1 x 10^^'' of PTB-CSF2 at optimum amplitude 
corresponds to 6P = 0.3 ppm. From figure HJ (5Pw is 7.7ppm so our measured bias 
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Figure 4. Calculated m = 1 DCP change in transition probability 5Pw of PTB- 
CSF2 as a function of microwave amplitude b for a single weakly coupled feed. 
The initial cloud offset parallel to the feeds is 1.7 mm. Solid line: tilt angle a\\ 
corresponding to maximum number of detected atoms. Dashed (dash-dotted) line: 
tilt angle — 2.54(+2.54) nirad with respect to the true vertical. 



suggests a suppression factor of > 25 and in turn that m = 1 DCP shifts are neghgible 
in figure [2l 

2.3. m = 2 DCP frequency shift 

The m = 2 quadrupolar cos(2(/)) phase variations are caused by transverse power flow 
from two feeds to the cavity walls. As for the m = 1 DCP shifts, the m = 2 DCP 
shift also vanishes for a horizontally symmetric cloud, launched vertically on the cavity 
axis and homogeneously detected. For a centered vertical launch, only the detection 
system can cause m = 2 DCP shifts [3]. But, if the detection laser beam propagates 
perpendicular to the optical axis of the imaging system and both are at 45° to the axis 
of the two cavity feeds, inhomogeneous detection will nonetheless be symmetric and 
produce no m = 2 DCP shift [3]. In PTB-CSF2 the orientation of the cavity feeds and 
the overall detection system are nearly optimal; the detection laser b 
37.5° from the feed axis, resulting in a suppression factor of 3.9 of the m = 2 DCP shifts 
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Figure 5. Calculated ra — 2 DCP change in transition probability 5P of PTB-CSF2 as 
a function of microwave amplitude b. Solid line: the initial cloud offset in the direction 
of the feeds is 1.7mm and the horizontal launch velocity is 0.7mm/s. Dashed (dash- 
dotted) line: the initial cloud offset in the direction of the feeds is 0.8(2.6) mm and the 
horizontal launch velocity is 0.7mm/s. Dotted line represents zero initial cloud offset 
and zero horizontal launch velocity (centered vertical launch). 

due to detection inhomogeneities. 

We evaluate the homogeneity of the detection system of PTB-CSF2 with a Monte- 
Carlo ray tracing simulation. We take into account the detection laser intensity 
profile and the geometric and optical parameters of the detection system. Because the 
detection laser intensity is high enough to reasonably saturate the fluorescence detection 
transition, the relative efficiency only varies from 100% to 90%. The fluorescence 
collection on the other hand varies much more. It has a minimum of 40% for atoms 
detected around the edges of the detection zone as compared to the fluorescence 
collection from atoms on the fountain axis. The dotted line in figure [5] shows the 
resulting calculated 5P for PTB-CSF2. 

The initial cloud offset also produces an m = 2 DCP shift. The total m = 2 DCP 
shift, including the detection inhomogeneity, is the solid line in figure |5l At optimum 
amplitude, h = 0.97, the m = 2 DCP frequency shift is 2.1 x 10^^^. The dashed and 
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Table 1. Individual DCP corrections m = 0, 1, 2 of PTB-CSF2 and their uncertainties 



(parts in 10^*^). 


Effect 


Correction 


Uncertainty 


DCP m = 


-0.018 


0.038 


DCP m = 1 (feed direction) 


-0.21 


1.08 


DCP m = 1 (perpendicular to feed direction) 




0.77 


DCP m = 2 


-0.21 


0.14 


total 


-0.44 


1.33 



dash-dotted lines represent transition probability changes for ±50% variations of the 
initial cloud offset. Considering the variations given by the dashed and dash-dotted 
lines as an upper limit for long-term drifts of the initial cloud offset, we get 1.4 x 10~^^ 
as the m = 2 DCP uncertainty. 

2.4- Summary of DCP frequency shifts and overall uncertainty budget of PTB-CSF2 

In table [1] the individual DCP contributions m = 0,1,2 of PTB-CSF2 are compiled 
and the overall correction is given, including its uncertainty. Our evaluation neglects 
DCP shifts from potential machining burrs [3l H]. Such shifts could be avoided by 
using a small cloud |S] or additional apertures to prevent atoms passing near the cutoff 
waveguide walls from being detected. 

We now include an entry in the uncertainty budget of PTB-CSF2 for the microwave 
lensing of the atomic wavefunctions by the microwave standing- wave O |9]. On the 
fountain ascent, the interaction of the atomic wave packets with the microwave standing- 
wave in the clock cavity gives rise to a spatially dependent phase shift that is nearly 
quadratic with transverse position, thus acting as a lens on the incident atomic wave [9] . 
The spatial phase variation from the dipole energy focuses one dressed state and 
defocuses the other during the first cavity passage. On the high frequency side of 
the central Ramsey fringe, the microwaves during the second cavity passage transfer the 
focused dressed-state to the final state of the clock transition and, on the low frequency 
side of the central fringe, the defocused dressed-state [9]. As a result, the apertures in 
the fountain cut more of the defocused dressed-state, producing an apparent positive 
frequency shift of the clock transition. 

The microwave lensing frequency shift is analytic if we ignore the apertures during 
the ascent, the microwave lens aberrations, and the detection inhomogeneities |5]. This 
analytic result predicts a shift of 10.9 x 10~^^. Including the lower aperture of the 
selection cavity during the ascent reduces the shift to 7.5 x 10"^^, because it truncates 
the large atomic cloud far from the cavity axis where the dipole forces are large. Adding 
the higher order lens aberrations gives a small correction, yielding 7.1 x 10^^''. Finally, 
with the detection inhomogeneity included [5], the predicted microwave lensing shift is 
8.3 X 10"^'' for PTB-CSF2. We take half of this correction as the uncertainty. 
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Table 2. Updated uncertainty budget of PTB-CSF2: Frequency shift type (origin), 
apphed frequency correction, and frequency correction uncertainty (parts in 10^^). 



Frequency shift 


Correction 


Uncertainty 


UaCli CLLIC ZjcCllldll hlilil 




yj.oo 


Blackbody radiation shift 


104.04 


U. 1 D 


Gravity + relativistic Doppler effect 


— O<3.0( 


U.UD 


Colhsional shift 


9.6 


3.0 


Cavity phase shift 


- 0.44 


1.33 


Microwave lensing 


- 0.83 


0.42 


AC Stark shift (hght shift) 


0.0 


0.01 


Majorana transitions 


0.0 


0.001 


Rabi puUing 


0.0 


0.002 


Ramsey pulHng 


0.0 


0.01 


Electronics 


0.0 


2.0 


Microwave leakage 


0.0 


1.0 


Background pressure 


0.0 


0.5 


total 


- 919.2 


4.1 



The updated uncertainty budget of PTB-CSF2 and the apphed corrections are given 
in table [2J For some corrections, i.e. quadratic Zeeman shift, blackbody radiation shift, 
colhsional shift, typical values are given. The uncertainty due to the colhsional shift is 
improved by a factor of ~ 2 compared to the first evaluation [7j by taking advantage of 
a better measurement statistics. 

The m = DCP shift accounts for our previously observed "microwave power 
dependence" so we no longer include this entry in the uncertainty budget of PTB-CSF2. 
We note that the currently largest uncertainty is due to collisions and is essentially given 
by its statistical measurement uncertainty, which will be reduced in future extended 
evaluations. Using rapid adiabatic passage [20] or tuning the collision energy to cancel 
the shift [21] could offer more profound reductions. 

3. Conclusions 

We theoretically and experimentaUy evaluate the distributed cavity phase shifts of the 
PTB-CSF2 fountain clock. The DCP frequency shifts due to m = longitudinal phase 
gradients are measured and compared with ah initio calculations. The experimental data 
agrees well with the model that shows a very small correction at optimum microwave 
amplitude. The m = 1 transverse phase gradients have been evaluated experimentally. 
At optimum microwave amplitude the gradient in the direction of the cavity feeds 
yields the largest DCP uncertainty contribution. The uncertainty is dominated by the 
statistical uncertainties of the measurements and the knowledge of the initial position 
and launch angle of the atomic cloud. The contribution from m = 2 transverse phase 
variations was modeled, accounting for the detection inhomogeneities and the initial 
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position and horizontal velocity of the atom cloud. We anticipate further reduction of 
the DCP uncertainty with better centering of the initial position of the atomic cloud. 
We evaluate the systematic frequency shift due to microwave lensing. Combined with 
a smaller uncertainty of the collision shift, the inaccuracy of PTB-CSF2 improves to 
4.1 X 10"^^ 
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